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Abstract 

We consider non-perturbative solutions of the Weinberg-Salam model at finite tem- 
perature. We employ an effective temperature-dependent potential yielding a first order 
phase transition. In the region of the phase transition, there exist two kinds of static, 
spherically symmetric solutions: sphalerons and bubbles. We analyze these solutions 
as functions of temperature. We consider the most general spherically symmetric fluc- 
tuations about the two solutions and construct the discrete modes in the region of the 
phase transition. Sphalerons and bubbles both possess a single unstable mode. We 
present simple approximation formulae for these levels. 
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1) Introduction 

The Weinberg- Salam model describes the electroweak interactions in the pertur- 
bative regime with remarkable accuracy. The model also predicts non-perturbative 
processes, such as instanton or sphaleron induced transitions, leading to baryon number 
violation. The suggestion [1] that the baryon asymmetry of the universe may possibly 
be explained within the standard model or its variants, recently created much interest 
in electroweak sphalerons and bubbles, non-perturbative solutions of the model. 

Sphalerons [2] are relevant for baryon number violating processes in the early uni- 
verse [3-5]. At a given temperature the rate of baryon number violation is largely 
determined by the three-dimensional action of the sphaleron solution, which enters in 
the Boltzmann factor, while the preexponential factor contains the frequencies of small 
oscillations about the sphaleron [3] . If the rate is large enough, the baryon asymmetry 
of the universe may be generated during the electroweak phase transition. 

However, for the creation of the baryon asymmetry of the universe during a first 
order electroweak phase transition, the details of this phase transition turn out to be 
important and need to be better understood [6-11]. The first order phase transition pro- 
ceeds via the creation of bubbles of the new broken phase embedded in the surrounding 
unbroken phase. The rate for the creation of bubbles is largely determined by the three- 
dimensional action of the bubble solution, which enters in the Boltzmann factor, while 
the preexponential factor contains the frequencies of small oscillations about the bubble 
[7]. 

Bubbles and sphalerons both arise as non-perturbative solutions of the finite tem- 
perature electroweak field theory. In the effective potential approach [12], the "Mexican 
hat" Higgs potential of the electroweak model is supplemented by temperature depen- 
dent terms. It is the presence of a term cubic in the Higgs field which renders the phase 
transition first order. 

The sphaleron of the electroweak theory, present at zero temperature, persists at 
finite temperature up to the region of the phase transition [13]. In contrast, the bubble 
solution exists only at finite temperature, in the region of the phase transition [7] . Here 
we study these non-perturbative solutions of the temperature dependent model in the 
region of the phase transition. In the first part of the paper we discuss the properties 
of these solutions, e. g. we address their three-dimensional action, their shape and 
their size. In the second part of the paper we analyze the discrete modes about these 
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solutions. We consider the most general spherically symmetric fluctuations, diagonalize 
the quadratic form in these fluctuations and compute the eigenvalues of the normal 
modes. We flnd that both sphalerons and bubbles possess a single unstable mode. For 
the negative eigenvalues we present approximation formulae. We relate the discrete 
modes of the bubble to those of the antikink. 



2) Lagrangian and effective potential 

We consider the electroweak theory in the limit of zero mixing angle. The bosonic 
part of the Lagrangian reads 

L = + {D^^YiD^^) - V{<p, T) , (1) 

(f){x) = ^($+$)^ , 
with the covariant derivative and fleld strength tensor 

D,^ = id,-i^r'^K)^ , (2) 

F;, = d^Al - d.A; + geatcA^Al . (3) 
T) denotes the temperature dependent effective potential [4,9-11] 

T) = - + It'^' - STcl>' + c , (4) 

1 - t:?2 ' y^i 



4y2 
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(6) 



V is the Higgs field expectation value at zero temperature, and Mw, Mz, Mf denote 
the masses of the W"^ bosons, the Z° boson and the top quark, and c is an adjustable 
constant. 

The effective potential (4) gives rise to a first order phase transition, characterized 
by three critical temperatures. For low values of the temperature T, the minimum 
of V{(j),T) is attained at some value {(j){T)) ^ 0, while = corresponds to a local 
maximum of V{(j), T). At T = Tc with 

r| = ^ (7) 

7 
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the second derivative of the potential at the origin changes sign. For T > Tc the 
extremum = turns into a local minimum separated from the absolute minimum 
{(j>{T)) by a small energy barrier. The two minima become degenerate at T = Tb with 



T, 
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For T > Tb the extremum = constitutes the absolute minimum of the potential. 
Finally at T = the local minimum disappears. 

To study temperature dependent quantities in the region of the phase transition it 
is convenient to represent the temperature via the variable ^ [10] 

Tc, Tfo, and Ta then correspond to ^ = 0, ^ = 2, and C = f > respectively. 

Throughout the paper we employ for the masses the values Mw = 80 GeV, Mz = 
92 GeV, Mt = 120 GeV and Mh = 45 GeV, and the vacuum expectation value for the 
Higgs field (0(0)) = v = 246 GeV, leading to the critical temperatures 

Tc = 0.284^; , T^ = 0.292v , T^ = 0.293^; . (10) 



3) Non-perturbative solutions 

In order to construct non-perturbative solutions of the effective lagrangian (1) we 
choose the gauge Aq — and we employ for the fields the static, spherically symmetric 
ansatz 

$(r) = -^L(r)(^j) , A^{r) = ^e,,an, (H) 

where L(r) and G{r) are radial functions. The three-dimensional effective action is then 
given by the functional 

S{T) = ^^^^ / dxx''a{x) , (12a) 
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with the action density 



the effective potential 
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and the dimensionless variable 

X = MwT , Mw = -zQV . 

The constant c in the effective potential has to be adjusted such that the integrand 
(12b) approaches zero asymptotically for the non-perturbative solution considered. 

Sphalerons 

The sphaleron, a well-known non-perturbative solution of the effective action (12) 
for T = 0, can be continued to finite temperatures, < T < T^. The sphaleron functions 
satisfy the boundary conditions 

G'(O) = 2 , G'(oo) = , (13a) 
L(0) = , L(oo) = -^^^ , (136) 

V 

where {(j){T)) denotes the non-trivial minimum of V{(j),T). The profiles of G{x) and 
L{x) are presented in Fig. 1 for T = Tc and T = 0. The action density (12b) of the 
sphaleron is presented for the temperatures T = T^, T = Tc and T = in Fig. 2. 
Both figures indicate that the sphaleron radius increases from r 3M^^ at T = to 
r ^ 9M^^ for T = Tfo. 

In Fig. 3 we exhibit the three-dimensional action of the sphaleron as a function of 
^ in the range of the phase transition. The three-dimensional action of the sphaleron 
follows an approximate scaling law, 

<5sp(T) = 53p(0)|||. (14) 

For an effective potential without the cubic term this scaling law is exact. 'S'gp(T) is 
well approximated by the scaling formula (14) [13]. The values predicted by the scaling 
formula exceed the exact numerical values only by a few percent, from 3% at T = Tc 
up to 6% at T = Tfo. For comparison the values obtained from the scaling formula (14) 
are also shown in Fig. 3. Since the scaling formula works so well, we have introduced a 
scaling factor 

im)) 

in Fig. 2 for the coordinate x and for the action density a{x). In the case of exact 
scaling, the three curves would be equal. 
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Bubbles 

In the range of the phase transition Tc < T < Tf), the equations of motion obtained 
from the functional (12) admit another non-perturbative solution, the bubble. It is 
through the formation of bubbles, that the first order electroweak phase transition 
proceeds. 

For the bubble solution, the gauge field consistently vanishes, G{x) = 0, and the 
equation for the Higgs field function L{x) reads 



dx^ X dx 

Choosing c = 0, a regular solution with finite three-dimensional action is obtained, 
where the function L{x) satisfies the boundary conditions 

dL 

— (0) = 0, L(oo)=0, (16) 

to be constrasted with boundary condition (13b). Note, that Derrick's theorem, forbid- 
ding the existence of static solutions of finite energy of a scalar theory in three spatial 
dimensions, applies only to a positive definite potential. 

The phase transition occurs in the temperature interval 1.3 < ^ < 1.8 [10,11]. The 
bubble function L{x) is presented in Fig. 4 for the temperatures ^ = 1.3, ^ = 1.8 and 
^ = 1.9. Only when ^ > 1.6 the bubble starts to develop an interior region, where L{x) 
is approximately constant. When ^ = 1.8 the size of this interior region is only about 
the size of the wall, and roughly 60M^^ wide. When ^ = 1.9 the interior region is only 
about twice the size of the wall. Increasing the temperature beyond ^ = 1.9 leads to a 
rapid increase in the size of the interior region, while the wall approaches its asymptotic 
shape, an antikink solution [10,11] 

LW^f [l-ta„h(-i=|^(.-X))]. (17) 

where X denotes the location of the antikink, which in the limit ^ = 2 moves towards 
infinity. 

Throughout the region 1.3 < ^ < 1.8 the thin wall approximation cannot be em- 
ployed. This is also clear from Fig. 5, where we show the action densities for the 
temperatures ^ = 1.3, ^ = 1.8 and ^ = 1.9. For the potential considered, the thin wall 
approximation is valid only for ^ very close to 2. 
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In Fig. 6 we show the quantity X, defined as the position where the derivative of 
the Higgs field function L'{x) assumes its maximal absolute value 

X : L\X) < , L"{X) = . (18) 

This value of X we employ as a definition for the location of the bubble wall. It is 
consistent with the use of X in eq. (17) in the limit ^ = 2. Approaching this limit we 
observe a rapid increase in X. Here X can be interpreted as the bubble radius. 

Fig. (6) also exhibits the quantity D 

which we interpret as a measure for the size of the wall. In the limit ^ 
approaches the corresponding analytical value of the antikink 

^ ^ 2MhMw 
STbV 

The value of the Higgs field function at the origin, L{0), is shown in Fig. 7 as a 
function of ^. L(0) starts from zero at ^ = 0, and increases with increasing ^ until about 
^ = 1.7, where it gets close to (0(T)). Then it decreases along with (0(r)), approaching 
(0(T)) from below. 

The three-dimensional action of the bubble obtained from eqs. (12) with c = is 
presented as a function of ^ in Fig. 3. We find good agreement with the semi-analytic 
formula presented by Adams [14]. In the limit ^ = 2, the action diverges, and the bubble 
ceases to exist. 

4) Normal modes 

Sphalerons and bubbles are extrema of the effective action (12). None of them 
is a stable minimum. To exhibit their instabilities one needs to consider general fluc- 
tuations about the extremal conflgurations and to construct the parts of the action 
functional quadratic in these fluctuations. The quadratic parts then represent operator 
equations, whose discrete negative eigenvalues correspond to the unstable modes of the 
non-perturbative solutions. 

By construction, the non-perturbative solutions discussed above are invariant un- 
der S0{3)d, a diagonal subgroup of SU{2)i x 50(3)^ (i. e. the group of weak isospin 



(19) 

2 D smoothly 
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and the group of spatial rotations). General fluctuations about the non-perturbative 
solutions can be decomposed according to the irreducible representations of S0{3)d, 
whose contributions decouple in the quadratic form. Here, we consider only the sector 
of spherically symmetric fluctuations, where the instabilities of the zero-temperature 
sphaleron are known to reside as well as the unstable breathing mode of the bubble. 

Normal modes about the sphalerons 

The most general spherically symmetric fluctuations about the sphaleron consists 
of flve radial functions 



A°;{x) = ^^^ (G{x) + ijGix)]eibaXb + i^H{x){5ia - XiXa) + V2^Kix)XiXa ■ (21 

ax L\ / J 



My 

gx 

The fluctuations if^ and i/jq decouple in the quadratic form from the fluctuations x, i^H 
and ifjK [15], yielding a system of two equations (2-channel) and a system of 3 equations 
(3-channel), respectively, 

/// fr2 3(G- 1)- - 1 2\ / V2LG , 

\ X / X ^22) 



V2LG , fG^ 4Arf2y 2X ^ 



1 , V2(G'- 1) , 

= ~-y^K — WH + Lx 



2V2{G-l-xG') , /2 3(^-1)2-1 2X V2LG 

\ X"^ / X 

For T = 0, the negative mode of the sphaleron, present in the 3-channel [15], is 
unique for Mh < l2Mw [16,17], while new directions of instability develop for Mh > 
12Mw The 2-channel possesses a single discrete positive mode for < Mh < l.SMv^ 
[17]. 

Also for T > the sphaleron possesses (at Mh = 45 GeV) a unique negative mode 
in the 3-channel, shown in Fig. 8. Analogous to the approximate scaling behaviour 
(14) of the three-dimensional action we find an approximate scaling behaviour of the 
eigenvalue a;^ of the negative mode 

<(r)=<..»i|^|^, (24) 
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which again is exact for an effective potential without cubic term. Using a;^(0) = —1.9, 
the deviation of the numerical values of a;^(T) from the approximate scaling values (24) 
is on the order of 10% at T = Tc (i. e. ^ = 0) and increases to about 20% at T = 
(i. e.e = 2). 

Normal modes about the bubbles 

For the discrete modes of the bubble the 2-channcl cqs. (22) is the relevant channel. 
For the bubble the two equations of the 2-channel decouple, since G{x) — 0. The 
equation for the fluctuation i/jq has a repulsive potential and therefore no normalizable 
discrete modes. 



The equation for the fluctuation il; reduces to 

9' 

The potential in eq. (25) is "deep enough" as to allow for a state with negative eigenvalue 
(i. e. a;^ < 0). This eigenmode represents the unstable breathing mode of the bubble, 
which lets bubbles smaller than the critical bubble shrink and larger bubbles grow. 
The negative eigenvalue a;^ is presented in Fig. 9. It starts from zero at ^ = 0, where 
the bubble solution appears, and ends at zero at ^ = 2, where the bubble solution 
approaches the antikink eq. (17), which also has a zero eigenvalue. For ^ — 2 the mode 
of the bubble indeed approaches the zero mode of the antikink 

Fig. 9 also shows approximate values for the eigenvalue a;^, obtained from the 
formula 

= , (26) 

where X is deflned in eq. (18). This approximate formula works amazingly well through- 
out the full region < ^ < 2, and not only in the vicinity of ^ = 2 [18], where X 
symbolizes the radius of the bubble. 

Beside the negative mode, eq. (25) allows for a discrete positive mode, also exhibited 
in Fig. (9). For small values of ^ the positive eigenvalue o;^ increases almost parallel to 
the continuum deflned by 

-)■ 



2 4A/T 



10 



When ^ approaches the limiting value of 2, however, uP' bends away from the continuum 
line. The positive mode of the bubble now approaches the positive mode of the antikink 



For the 3-channel of the bubble the equations do not decouple. Our analysis indi- 
cates, however, that there are no discrete modes in this channel. 

5) Conclusion 

In this paper we have studied the sphaleron and the bubble, two non-perturbative 
solutions of the electroweak model in the region of the electroweak phase transition. The 
temperature dependence of the model is taken into account via an effective potential, 
which assures a first order phase transition. 

The sphaleron solution, its action density and in particular its three-dimensional 
action follow an approximate scaling behaviour with respect to the well-known zero 
temperature sphaleron. The action decreases with the scaling factor A = (0(T))/(0(O)), 
while the size increases with this scaling factor. At the temperature the sphaleron is 
bigger by a factor of about 2 and at T5 it is bigger by a factor of about 3 than the zero 
temperature sphaleron, giving it a size of about 10 — 20Myy in the region of the phase 
transition. 

The sphaleron has a single unstable mode in the region of the phase transition. 
The negative eigenvalue also shows an approximate scaling behaviour. Indeed, the 
eigenvalue o;^ of the negative mode at finite temperature is rather well approximated 
by the negative eigenvalue at zero temperature scaled by the factor A^. 

The bubble solution only exists for sufficiently large temperatures, Tc < T < Tb. 
In the limit T ^ the solution converges to the trivial solution (i. e. all fields are 
vanishing). In the limit T — > T5 the bubble solution approaches the antikink solution, 
where the bubble wall moves towards infinity, and the three-dimensional action diverges. 

In the vicinity of the phase transition the bubble function develops an interior 
region where it is approximately constant, followed by the region of the bubble wall. At 
temperatures where the phase transition occurs, the size of both regions is of the same 
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order of magnitude. The thin wall approximation is not justified here. It is justified 
only for T — > T^, where the size of the interior region rapidly increases, while the wall 
approaches an asymptotic shape and size, given by the shape and size of the antikink. 
The size of the bubble wall, being a few times the size of the sphaleron, is thus of the 
same order of magnitude as the size of the sphalerons at these temperatures. 

The bubble has a single negative mode, the breathing mode. The frequency of this 
mode starts from zero at Tc and again approaches zero at T^. For T — > the negative 
mode of the bubble and its eigenvalue approach the zero mode of the antikink. The 
negative eigenvalue is very well approximated by the simple formula uP = —2/X^, where 
X is defined as the value, where the bubble function changes most strongly, i. e. |L'(X)| 
is a maximum. For large temperatures, X may be viewed as the bubble radius. 

The bubble also possesses a positive discrete mode. As a function of temperature 
the eigenvalue of this mode increases parallel to the continuum. However, for T ^ the 
eigenvalue bends away from the continuum: the positive mode of the bubble approaches 
the positive mode of the antikink. 
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7) Figure captions 

Fig. 1 The gauge field function G{x) and the Higgs field function L{x) of the sphaleron are 
shown as a function of the dimensionless coordinate x = Mwf for T = (dashed 
curves) and T = Tc (solid curves) . 

Fig. 2 The action density (7{x) of the sphaleron is shown as a function of the coordinate 
Ax for the temperatures T = (dashed), T = Tc (dotted) and T = (solid). A is 
the scaling factor {(f){T)) / {(f){Tb)) . 

Fig. 3 The action S{T) of the sphaleron and of the bubble in units of Mw/ctw are shown 
as functions of temperature for Tc < T < Tb, expressed in terms of the variable 
^ for < ^ < 2 (solid). For comparison the approximate scaling values for the 
sphaleron action are also shown (dashed). 

Fig. 4 The Higgs field function L{x) is shown for the bubble as a function of the dimen- 
sionless coordinate x = Miyr for the temperatures ^ = 1.3 (dotted), ^ = 1.8 (solid) 
and ^ = 1.9 (dashed). 

Fig. 5 The action density a{x) is shown for the bubble as a function of the dimensionless 
coordinate x = M^r for the temperatures ^ = 1.3 (dotted), ^ = 1.8 (sohd) and 
^ = 1.9 (dashed). 

Fig. 6 The coordinate X, where the bubble derivative L'{x) has its maximal absolute 
value, and the quantity D = —L{0)/L'{X) are shown as functions of the tempera- 
ture expressed in terms of ^. 

Fig. 7 The quantity L(0) is shown for the bubble as a function of the temperature ex- 
pressed in terms of ^ and compared with the minimum of the potential {(f){T))/v. 

Fig. 8 The eigenvalue of the negative mode of the sphaleron is shown as a function of 
the temperature expressed in terms of ^ (solid). For comparison the approximate 
scaling values for the eigenvalue are also shown (dashed). 

Fig. 9 The eigenvalues of the discrete modes of the bubble are shown as functions of 
the temperature expressed in terms of ^ (solid). For comparison the approximate 
formula o;^ = —2/X'^ (dashed), and the continuum (dotted) are also shown. 
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